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Abstract. We explore the possibility to measure the elastic and inelastic uiN cross section in p+d^d+u+p sp and 
p+A reactions. Our studies indicate that the elastic scattering cross sections can be determined for u> momenta above 
1 GeV/c in p+d reactions by gating on high proton spectator momenta whereas the uiN absorption cross section down 
to low relative ui momenta is most effectively studied in p+A reactions at beam energies 2.0-2.7 GeV. 

PACS. 13.75.Cs Nucleon-nucleon interactions - 25. 40.-h Nucleon-induced reactions - 25. lO.+s Nuclear reactions 
involving few-nucleon systems - 11.80.Fv Eikonal approximation - 24.10.Lx Monte Carlo simulations 



1 Introduction 

The w-meson properties at finite nuclear density have become 
of recent interest especially in the context of dilepton studies [jll 
gj, where one hopes to see a shift of the meson pole e.g. at den- 
sity po~0.16 fm -3 []3j^,g,g[j. However, the in-medium prop- 
erties of the vector mesons - as reflected in their spectral func- 
tions - are presently a matter of strong debate. Whereas at high 
baryon density and/or temperature some scaling of the meson 
poles with the scalar quark condensate <qq> might be ade- 
quate [0], the properties of the mesons in a hadron gas at low 
density are essentially determined by the scattering amplitudes 
with the hadrons in the local environment [bJ. By means of 
dispersion relations the real and imaginary part of these scat- 
tering amplitudes - in the low density approximation - directly 
convert to mean-field potentials and scattering rates, respec- 
tively [||] . Thus it is of profound interest to obtain more precise 
information especially about the loN cross sections to construct 
its spectral function at low baryon density. 

We recall that the strong coupling of the p-meson to nu- 
cleon resonances leads to large inelastic cross sections with nu- 
cleons at low relative momenta which in turn results in a large 
collisional width at density p a [|,|,0[n],|iJ[l3|]. If this melt- 
ing of the p-meson will also occur for the u> meson is of ma- 
jor importance for dilepton experiments, that aim at detecting 
moderate but clear peaks in the dilepton invariant mass spectra. 
Though there are a couple of theoretical estimates for the uoN 
cross sections [0[T|op],pl[i]| the theoretical predictions 
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widely differ such that related experiments will have to clarify 
the situation. 

Moreover, the amplitude for w-meson production in proton- 
nucleon reactions close to threshold so far is scarcely known - 
especially with a deuteron in the final state - and of interest in 
itself. 

In this work we aim at contributing to the latter question by 
investigating theoretically within the Multiple Scattering Monte 
Carlo (MSMC) approach and transport model calculations, to 
what extent we might obtain experimental information from 
various reactions at existing accelerators that can be performed 
in the near future. In Section 2 we will study in detail the per- 
spectives of the reaction p + d^d+uj+p sp close to threshold. 
In Section 3 proton induced reactions on nuclei will be inves- 
tigated within a coupled channel transport approach to explore 
the production and propagation of the w-meson in finite nuclei 
whereas Section 4 is devoted to a summary of our work. 



2 Effects from u> rescattering in the reaction 

p+d^d+UJ+p sp 

In Ref. [ |l7| ] we have considered the cross section of the reaction 
p+d-^d+uj+p within the framework of the spectator mecha- 
nism close to threshold energies (see diagram a) in Fig. pi). This 
mechanism gives the dominant contribution to the cross section 
if the momentum of the proton-spectator from the deuteron is 
below 100-150 MeV/c. The amplitude corresponding to the di- 
agram a) of Fig. [l] can be written in the following form in the 
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deuteron rest frame 



f(pn -» dV) / d 3 r exp(-ip 2 r) (j> d (r), (1) 



M a 



where V is the vector meson (u>), P2 is the proton spectator mo- 
mentum, </>d(r) is the deuteron wave function and f(pn— *dV) 
is the amplitude of the reaction pn-^dV, related to the differ- 
ential lab. cross section as 



da(pn — > dV) 
dM 



= \f(pn->dV)\ 



n-". 



V 



a) 




(2) 



b) 



Fig. 1. The diagrams for vector meson (V) production in the 

p+d— *d+V +p sp reaction without (a) and with VN rescattering (b). 



The probability for the produced meson to rescatter can be 



found by integrating the term proportional to 
equal to 



I Mb I 2 , which is 



W = a e i{Vp->Vp) dzG(z)\4> d (b = 0,z)\' 



(6) 



Now the momentum distribution of the spectators can be 
written in the form 



dN(p) 
p 2 dp 



=N(\ip d {p)\ 2 + \^ resc {p)\ 2 +InterferenceTerm), 



(7) 

where the first term is the contribution of the spectator mecha- 
nism, with ipd(p) denoting the deuteron wave function in mo- 
mentum space, while N is the normalization factor. The sec- 
ond term |tfV e sc(p)| 2 appears due to the rescattering of all par- 
ticles in the final state on the proton-spectator added incoher- 
ently with the probabilities defined by Eq. (Q). The interference 
between the spectator and rescattering amplitudes is only im- 
portant in a narrow region of spectator momenta, where both 
contributions are of the same order of magnitude. For a more 
detailed discussion of this point see, e.g. Refs. [ |l8[[l9| , po[pi[ |). 
In order to obtain closer information on the vector meson 
rescattering we employ the Multiple Scattering Monte Carlo 
(MSMC) approach. An earlier version of this approach, de- 
noted as Intra-Nuclear Cascade (INC) model, has been applied 
to the analysis of n and oj production in pA and pA interac- 
tions in Refs. [E2LG3I], respectively. Recently this version of the 
INC model has been extended to incorporate in-medium modi- 
fications of the mesons produced [Q,Q[| in hadron-nuclear colli- 
sions. For an explicit presentation of the production cross sec- 
tions employed we refer the reader to Ref . [B[| or the review pi] . 



The produced vector meson can rescatter elastically on the 
proton-spectator, thus transferring momenta larger than 200 
MeV/c. Near threshold the momentum of the produced oj me- 
son is quite large: about 0.8-1 GeV/c. Therefore its wave length 
is much smaller than the average distance between two nucle- 
ons in a deuteron. In this case we can describe the rescattering 
amplitude (see diagram b) in Fig. [y) within the framework of 
the eikonal approximation (cf . Ref. [[18]] ) 



M b = -f(pn -► dV) J d 3 r ex P Hp 2 r) 0(z)r(b)M^), 

(3) 
where the z- axis is directed along the V-meson momentum. 
The elastic Vp - scattering amplitude is related to the profile 
function r(h) by the standard expression 



f(V P -> Vp) 



1 



2nik 



d 2 b exp(-iqb) T(b) (4) 



with k denoting the lab. momentum of the V-meson and q the 
momentum transfer. 

The probability to detect the deuteron at the solid angle 
dfid in coincidence with the spectator of momentum p 2 is re- 
lated to the differential cross section as 



d 5 a(pd —> dVp 2 ) 
df?d 3 p2 



\M a + M b 



(5) 
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Fig. 2. Momentum distribution of proton-spectators in the reaction 
pd— »37r + 27r~p. The data are taken from Ref. ||2a] while the solid line 
is calculated using Eq. (bl). 
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However, the INC model is valid only for medium and heavy 
nuclei and cannot be directly applied to deuterons. In order to 
perform simulations of scattering events in the case of hadron 
- deuteron interactions we use the Monte Carlo approach for 
the description of single and double scattering terms. The mo- 
mentum distribution of a nucleon in the deuteron is described 
by the Fourier transformed deuteron wave function squared 
\^pd(p)\ 2 - The total probability of rescattering is calculated for 
each secondary particle using Eq.(6). All events are generated 
using realistic angular distributions for the elementary subpro- 
cesses. In the case of u> meson rescattering this information is 
not available from data and we employ different assumptions 
on its cross section, i.e. in magnitude and angular distribution. 

Experimentally, proton-spectator spectra have been mea- 
sured in pd annihilation p4| , 25J] to some extent. As an exam- 
ple - and to testify the applicability of the model to such type 
of reactions - we present in Fig. || the calculated momentum 
distribution of proton spectators in the reaction pd— >37r + 27r~p 
at rest in comparison to the data from Ref. [|25||. These data 
are easy to interpret: the sharp peak with a maximum near 60 
MeV/c corresponds to the spectator mechanism and the long 
tail above 300 MeV/c is due to the rescattering of pions. The 
solid line in Fig. || is calculated using Eq. (Q) (without the in- 
terference term), where the deuteron wave function was taken 
in the parameterization of Ref. Pq ] with S- and D-waves in- 
cluded. The dependence of the second term on the momentum 
p in the lab. system is determined by the kinematical conditions 
(like the mass of the rescatterred particle) and by the angular 
dependence of the pion-nucleon scattering amplitude. We find 
that the spectrum calculated by the MSMC approach is in good 
agreement with the experimental data. 

The same approach was used to simulate the momentum 
distributions of spectator protons in the reaction p+d^d+iv + 
p sp . Since the elastic rescattering cross section of the w-meson 
is basically unknown and especially its angular distribution at 
high relative momenta, we study the observable effects within 
different models. In Fig. |] we present three assumptions for 
the angular distribution of u)N elastic rescattering (upper left). 
These assumptions correspond to a parameterization of the dif- 
ferential cross section da/dt — Cexp(6i) with slope parame- 
ters b = 0, 6 and 12 GeV 2 , respectively. 

The resulting spectator momentum distributions are shown 
in the upper right part of Fig. g. The lowest (dotted) histogram 
describes the contribution of the spectator term when rescat- 
tering is absent. In this case the proton momentum spectrum 
only extends to 450 MeV/c. The most pronounced effect from 
rescatterings is seen for the isotropic model (solid histogram) 
where the proton momentum spectrum extends to 1 GeV/c. 

Furthermore, the distributions in the azimuthal angle <P be- 
tween the scattering planes p—p^ and p—p sp for spectators 
with p sp <150 MeV/c (middle left), 150<p s)9 <300 MeV/c 
(middle right), 300<p sp < 450 MeV/c (lower left) and for fast 
spectators with p sp > 450 MeV/c (lower right), show the most 
pronounced rescattering effects for fast spectators in line with 
the upper right part of the figure. 

It is clearly seen that the spectator mechanism dominates 
at p sp <30Q MeV/c while at p sp >450 MeV/c the main con- 
tribution stems from the rescattering term. The best way to 
select the rescattering term is to look for the correlation be- 



tween the azimuthal angles of the two scattering planes p—p^ 
and p—psp (lower figures). This correlation is much more pro- 
nounced for fast spectators (p sp >450 MeV/c) than for slower 
ones (p sp <300 MeV/c). It is worth to point out that the slope 
of the spectator spectrum depends on the slope of the lo— p an- 
gular distribution, while the correlation in the angle <P depends 
mainly on the value of the u> N elastic cross section. 
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Fig. 3. Angular distribution of u)N elastic rescattering (upper left) 
within 3 different models indicated by the solid, dashed and dotted 
histograms (see text). The resulting proton-spectator momentum dis- 
tributions are shown in the upper right part for the 3 models in the 
reaction p+d— >d+u)+p sp . In all calculations the loN rescattering 
cross section u e ;=40 mb was assumed. The distributions in the az- 
imuthal angle $ between the scattering planes p—pu, and p—p 3p for 
spectators withp sp <150 MeV/c (middle left), 150<p sp <300 MeV/c 
(middle right), 300<p sp <450 (lower left) and for fast spectators with 
p 3p >450 MeV/c (lower right) are shown additionally. 
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Thus performing a cut for p sp >250 MeV/c and 150° < 
<&(u>, p sp )<210° one can determine essentially the overall mag- 
nitude of the ujN elastic cross section (cf. upper part of Fig. Q), 
whereas the number of events as a function of the cut in the 
proton-spectator momentum p sp provides information on the 
qj—N angular distribution. This is demonstrated in the lower 
part of Fig. [i] for the 3 models of the uj— p angular distribu- 
tions, where the solid, dashed and dotted lines correspond to 
the upper left part of Fig. [3j 



150°<*(w,p sp )<210° 




Fig. 4. The number of events for the cut ps P >250 MeV/c and 
150°<(/)(cJ,p S p)<210 o as a function of the elastic cup cross section 
for 3 different angular distributions (solid, dashed and dotted lines) 
from the upper left part of Fig. [3j is displayed in the upper part of the 
figure. The lower part shows the number of events as a function of the 
cut in the proton-spectator momentum p 3p for the 3 different angular 
distributions assuming an overall magnitude of a e i (up) =40 mb in the 
reaction p+d^d+ui+p sp - The corresponding slope parameters b are 
given in GeV~ 2 (see text). 



It is worth to point out that the relative contribution of the 
rescattering term is about a few percent as compared to the 
spectator term. In the case of small spectator momenta, when 
the spectator term is dominant, the reaction p + d^d+uu +p sp 
can be identified by the missing mass method measuring the 
forward deuteron in coincidence with a slow spectator proton. 
We note that the spectator detection with thin semiconductor 
counters has already successfully been applied for ?/-meson 
production in pd collisions at CELSIUS p7|,|28|]. However, the 
selection of the reaction p + d^d+cu+p sp with fast (p sp >250 
MeV/c) rescattered spectator protons is more difficult since 
background particles (e.g. 7r-mesons) can be misidentified as 
spectator protons. Thus, the direct detection of the produced 
w-mesons by the detection of its decay products like 7r°7 is 
mandatory. We suggest to use a neutral particle detector and to 
gate on three energetic photons, where two photons should be 



correlated in their invariant mass to the ir° whereas the third 
photon should add to the invariant mass of the w-meson ac- 
cording to its Dalitz decay uj^it q +j with a branching ratio of 
-8.5%. 

The momentum and angular distributions of 7r°-mesons and 
photons in the reaction p+d^d+uj(— >7r 7)+p sp at T p =2 GeV 
are shown in Fig. pi The momentum spectra of pions and pho- 
tons have maxima at 0.5 GeV/c. Thus the detection of such 
energetic photons will be a very good trigger for u> production, 
i.e. the selection of the reaction can be done by measuring the u> 
with a fast forward deuteron in coincidence. We note that such 
measurements might well be performed at COSY-Julich. 
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Fig. 5. The momentum and angular distributions of n° mesons and 
photons from the reaction p+d—*d+u>(—yKj)+p 3 p at T p =2 GeV. 



3 Proton induced reactions on nuclei 

The production of vector mesons in proton induced reactions 
on light and heavy nuclei has been studied before for p-meson 
propagation in the medium by looking at 7r + 7r— invariant mass 
spectra in Ref . [ |29| ] . For the dynamical studies we again employ 
the coupled channel transport model [ p0| ] with elementary pro- 
duction cross sections as described in Refs. [pi. The transport 



approach [ 30 1 allows to account for the final state interactions 
(FSI) of all hadrons and especially the elastic and inelastic u>N 
scatterings. For the experimental detection of the cj-meson we 
again suggest to gate on three energetic photons. 
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In order to include the effects from uN elastic rescattering The momentum distribution of cj-mesons in beam direc- 

we employ the cross section from Ref. [ [lq ] in the parameteri- tion, 

zation „ 

da f 2 " a 

o- rf = 5.4 + 10 *exp(-0.6p ttf )[mb], (8) ~M\ ~ J P± ~ ^V 



where p u is the ujN relative laboratory momentum in GeV/c. 
For the inelastic channels we have to consider various reactions 
(cf. Ref. [[lq]). Among them the u + N^tt+N cross section 
can be determined quite well from the available data on the 
exclusive u> production in nN reactions from Ref. [ pl| ] by ex- 
ploiting detailed balance. The uiN inelastic cross sections due 
to other reaction channels have been calculated in Ref. Jlq | 
by a meson-exchange model and, in principle, are determined 
within the uncertainties of the coupling constants and the form 
factors employed. As a guide we thus use the inelastic ujN 
cross section <Ji ne i from Ref. Jlq ] within the parameterization 



(Jinei =20 + 4.0/p w [mb], 



(9) 



is shown in Fig. morp+ 12 C collisions at beam energies of 2.0, 
2.4 and 2.7 GeV. The dashed lines in Fig. ^ (full dots) show the 
calculations by taking into account the elastic and inelastic uoN 
interactions, while the solid lines indicate the results obtained 
without FSI (full squares). 

At all energies considered the lj longitudinal momentum 
distribution in the laboratory is extending over a broad momen- 
tum regime from negative momenta even up to pu ss 2.0 GeV/c 
at T; a f,=2.7 GeV. By comparing the dashed and solid lines we 
find a substantial reduction of the w-meson yield at medium 
and high momenta even for 12 C. Furthermore, for u momenta 
in the Fermi-momentum regime pf< 0.23 GeV/c a suppres- 
sion is less obvious. This is due to elastic ojN^uoN scattering 
which essentially populates low pu in the laboratory system. 



with p u given in GeV/c. 



p+ 12 C^w+X 



■ without FSI 
• with FSI 



T=2.0 GeV 




1.5 2 

Pu (GeV/c 

Fig. 6. The longitudinal momentum distribution of oj-mesons pro- 
duced inp + 12 C collisions at 2.0, 2.4 and 2.7 GeV. The symbols show 
the transport calculations while the lines are drawn to guide the eye. 
The squares show the calculations without ujN final state interactions 
while the full dots include FSI. 



p + ,2 C^w + X at 2.4 GeV 




transport model 
Glauber model 



p. (GeV/c) 

Fig. 7. The ratio « of the a>meson spectrum without FSI to the ui- 
spectrum including FSI (triangles) from the transport calculation for 
p + 12 C collisions at 2.4 GeV. The solid line is drawn to guide the 
eye while the dotted line shows a Glauber calculation including the 
momentum dependent luN cross section. 



The latter phenomenon is more clearly demonstrated in Fig. Uj 
where we display the ratio k of the u spectrum without FSI to 
the u) spectrum including FSI (triangles) from the transport cal- 
culation for p+ 12 C at 2.4 GeV. The solid line in Fig. [7] is drawn 
to guide the eye. Moreov er, since the G lauber approximation is 
traditionally considered 



33 34j,B3] as an effective way to 



evaluate the cross section of an unstable particle with nucleons 
we perform a comparison between the Glauber model and the 
results from the transport approach. 

We recall that within the Glauber model the A-dependence 
for w-meson production and propagation in the nucleus is given 
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+ 00 +OO 2-7T 

A e ff — / bdb / p(b, z)dz / d(f> 
-00 

z 

x[exp(-o- pA r p(b,£)d£-a uN f p(r c )d()], (11) 



where p(r) is the single-particle density distribution normal- 
ized to the target mass number. The last integration in Eq. ( |ll| ) 
proceeds over the path of the produced w-meson 

r C = (b + (cos(f)sin8) 2 + ((sin(f)sin9) 2 + (z + (cos9) 2 , (12) 

where 9 is the emission angle of the w-meson relative to the 
beam direction. Moreover, a P N is the pN cross section while 
cr^N is the total coN cross section. 

Within the Glauber model the pii -dependence of the factor 
K is given by the ratio of A e ff - calculated with the momentum 
dependent u)N cross section (Q) - to the result for <t w jv=0. The 
dashed line in Fig. [7] shows the factor k calculated by Eq. (pT[). 
It is clear that the calculations by the Glauber model fail in de- 
scribing the FSI as generated by the transport calculation since 
the Glauber analysis discards elastic rescattering and accord- 
ingly misses the population of the low momentum part of the w 
spectrum. Thus the factor k calculated by the transport model 
substantially differs from the Glauber calculations at low uo- 
meson momenta. Furthermore the difference at p» >500 MeV/c 
is due to the fact that the transport calculations account for the 
uj production by secondary production processes, while the cal- 
culations by the Glauber model do not. 

Indeed, we note that at subthreshold energies in p+A reac- 
tions the secondary production channel with an intermediate 
pion (p+N—tN+N+w, tt+N^uj+N) is the dominant one 
according to the studies performed in Refs. [g6,37|. On the 
other hand, the calculations by Eq. ([FlJ) account only for the 
direct pN-^pNuo production mechanism. Thus the application 
of the Glauber model is not appropriate in this case for an ex- 
perimental analysis. 

The mass dependence of the total production cross section 
of mesons can be exploited to provide some information on the 
strength of the meson absorption cross section. In this respect 
we show in Fig. |8J (full dots) the expected mass A-dependence 
of ui mesons from p+A reactions at beam energies T p =2.0, 2.4, 
and 2.7 GeV as calculated by the transport model without in- 
cluding FSI. The cross sections can be well described at fixed 
energy by a power law 



a(p + A -> uj + X) 



A a , 



(13) 



where the parameter a sheds some light on the reaction mecha- 
nism. The solid lines in Fig. Hshow a fit to the transport model 
calculations with parameters a indicated in the figure. Further- 
more, for the direct p+N^ui+p+N production (neglecting 
the w-meson FSI) we expect a«0.75 from the Glauber model 
(dashed lines in Fig. ph which, however, underestimates the 
mass dependence obtained from the transport model calcula- 
tions at all energies and especially at T p =2 GeV. 

As discussed in Ref. [po[]-[|3lfl the higher power in A is 
due to secondary reaction channels involving an intermediate 
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Fig. 8. The A-dependence of u>-meson production from p+A reac- 
tions at beam energy T p =2.0, 2.4, and 2.7 GeV calculated without 
FSI. The full dots show the transport model results while the solid 
lines are a fit by Eq. < |1 3f > with parameters a as indicated in the figures. 
The dashed lines show the calculations by the Glauber model without 
FSI, which provide a~0.75. 

pion, which are more frequent in heavy nuclear targets and at 
lower energy. On the other hand, at bombarding energies of 
2.4 and 2.7 GeV the primary reaction channel dominates such 
that we come closer to the results from the Glauber calculation. 
Furthermore, the difference between the parameter a obtained 
from the transport calculations to that from the Glauber model 
at beam energies T p =2.4 and 2.7 GeV actually indicates the 
contribution from the secondary reaction mechanism. 

The mass dependence of the oj production cross section in- 
cluding the FSI as resulting from the transport calculations is 
shown in Fig. ^ in comparison to the Glauber model employ- 
ing a constant luN cross section of 35 mb in the latter (dashed 
fines). The respective power law parameters a from the trans- 
port calculation differ with bombarding energy and approach 
a=2/3 above ~2.4 GeV bombarding energy characterizing the 
surface dominance of uj meson production in p+A reactions 
due to strong absorption. On the other hand, the Glauber model 
gives a rough description of the mass dependence at higher en- 
ergies for <7 w jv=35 m b (with a « 0.55). This sensitivity thus 
might also be explored experimentally, however, keeping in 
mind that the measurements should be performed at different 
bombarding energies to avoid systematic errors associated with 
different reaction mechanisms. 
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Fig. 9. The A-dependence of aj-meson production from p+A reac- 
tions at beam energy T p =2.0, 2.4, and 2.7 GeV calculated with u>N 
FSI. The full dots show the transport model results while the solid 
lines are a fit by Eq. (O) with parameters a as indicated in the fig- 
ures. The dashed lines show the calculations by the Glauber model for 
ct„jv=35 mb, which give a~0.55. 



4 Summary 

In this work we have explored the possibility to measure the 
elastic luN cross section in p+d^d+uj+p sp reactions and the 
inelastic ui N cross section in p+A collisions. Our studies based 
on MSMC and transport calculations indicate that the elastic 
scattering cross sections can be determined for lo momenta 
above 1 GeV/c in p + d reactions by gating on high proton 
spectator momenta whereas the uN absorption cross section 
is most effectively studied in p+A reactions around Ti a i,=2A 
GeV via the mass dependence of the total production cross 
section. Though the apparent power law in the target mass A 
is sensitive to the actual strength of the u N inelastic cross sec- 
tion, its actual momentum dependence cannot be extracted eas- 
ily from experimental data due to strong effects from elastic 
uN scattering processes. 
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